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ABSTRACT 

This  paper  presents  a  new  numerical  method  for  solving 
Schroedinger 's  equation  in  radial  form  in  order  to  get  a  set 
of  wave  functions  for  use  as  starting  values  in  the  solution 
of  the  Fock  equation.  As  a  necessary  adjunct  to  the  solution, 
an  improved  numerical  method  is  given  for  calculating  potentials 
for  use  in  the  Schroedinger  equation.   Finally,  a  method  is 
outlined  for  solving  the  Fock  equation  numerically,  using  the 
wave  functions  obtained  above  as  starting  values. 
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ENERGY  OF  THE  MOST  LOOSELY  BOUND  ELECTRON 
FOR  IONS  OF  COPPER,  SILVER,  AND  GOLD 

If  the  energy  and  distance  from  the  nucleus  of  the  most 
loosely  bound  electron  in  an  ion  of  a  metallic  lattice  can  be 
found,  it  is  possible  to  use  this  data  to  calculate  certain 
physical  constants  of  the  metal,  e.g.  the  elastic  constants. 

Wigner  and  Seitz1  have  shown  how  to  find  this  binding 
energy,  neglecting  exchange  forces,  by  solving  Schroedinger 's 
equation  in  radial  form: 


(D      -ii  <£*  .-[vM-eJe  -  o 


"2.  rw\    <A.    /" 

for  various  values  of  E.   It  is  found  that  for  each  such  value 
of  E  there  are  several  radii  at  which  the  condition 


(2) 


df  f 


is  satisfied.  Most  of  these  points  are  very  close  to  the 
nucleus  and  are  essentially  the  same  for  all  values  of  E„   If 
the  outer  two  radii  are  plotted  against  E  a  curve  results 
which  exhibits  a  minimum.   It  is  this  minimum  which  gives  the 

required  data. 

p 
Fuchs  has  extended  the  procedure  to  take  into  account 

exchange  forces.   His  method  involves  solving  the  Foci  equation 

1E.  Wigner  and  Seitz,  P.,  Phys.  Rev.,  4^,  804  (1955) 
2K,  Fuchs,  Proc.  Roy.  Soc,  A,  1£1,  585  (1955) 
^E.  Fock,  Zeits.  F.  Physik,  81,  195  (1955) 


.5 


(5)   -^1  <££  +[v(r)-6]n- J6(r,^^r')J.r'.  0  . 

The  most  practical  method  for  solving  this  equation  is 
first  to  solve  equation  (1)  and  then  use  the  wave  functions 
obtained  to  evaluate  the  integral  in  (5).  This  makes  the  right 
side  a  function  of  the  independent  variable  instead  of  zero. 
When  this  resulting  equation  is  solved  the  new  wave  functions 
are  then  used  to  re-evaluate  the  integral  and  the  procedure  is 
repeated  as  required. 

As  the  first  approximation,  then,  a  solution  of  (1)  is 
needed.  To  get  it,  it  is  necessary  first  to  calculate  the 

values  of  V  to  be  used. 

4 
From  the  work  of  Hartree  and  others  ,  tables  of  values  of 

Z,  the  effective  charge,  are  available  as  functions  of  /O  ,  the 

radius  measured  from  the  surface  of  the  nucleus.   These  values 

have  been  accepted  in  this  paper,  although  there  may  be  errors 

in  them  since  most  were  calculated  by  hand  before  the  era  of 

the  electronic  computer.   The  required  values  of  v,  however, 

were  recalculated,  where 

and  Z  is  the  effective  nuclear  charge  for  potential  which  is 
different  from  Z.  The  values  of  v  were  calculated  by 

%.   Hartree,  Proc.  Roy.  Soc,  A,  l4l,  282  (1955) 
A.  Douglas,  Hartree  D.,  and  Runciman,  W.  Proc.  Cambridge 

Phil.  Soc,  £1,  486  (1955) 
M.  Black,  Manchester  Memoirs,  ££,  29  (1954) 


integrating  the  field  intensity 


/ 


The  usual  method-?  of  calculating  v  is  to  solve  the 
differential  equation 

(6)  ^£  ..  gpl* 

<*e         e 

by  Euler's  method,  which  is  simple  but  inaccurate.   It  can  be 
shown  that  the  error  is  of  the  order  of 

(7)  6  '-    Ce'e^   (^% 


Some  simple  algebra  and  substitution  reduces  this  expression  to 

(8)  t-LJ^L\±^)-H?,)\ 

*p, 

since  Z  is  assumed  to  be  linear  between  any  two  adjacent  radii. 
Graph  1  gives  the  relative  %   error  for  Cu+  resulting  from  this 
calculation.   It  is  large  enough  so  that  a  more  sophisticated 
method  of  solution,  equation  (5) >  was  decided  upon.  The 
evaluation  was  done  using  a  Gaussian  5-P°int  quadrature  formula, 
which  is  based  on  the  assumption  that  between  any  two  neighboring 
radii  Z  is  approximated  by  a  third  degree  polynomial. 
To  evaluate  the  integral  in  (5) » 


(9) 

r 

5d.  Hartree,  Proc.  Cambridge  Phil.  Soc,  24,  114  (1928) 

6Kunz,  K.,  Numerical  Analysis,  McGraw-Hill,  N.Y. ,  1957, pgs  177-8 


it  is  first  necessary  to  normalize  the  limits  to  the  range  0  to  1 
in  order  to  use  readily  available  values  of  the  Gaussian  coef- 
ficients.  This  can  be  done  by  the  substitution 

do)  w.-  -zy-  -^rzrpT     +  '  I 

Equation  (9)  then  becomes 

(11)   A[f0=  4(f„-f,))e  u^Ofe-P,)"!^^"^ 

The  coefficients  for  the  Gauss  5-P°int  formula  are 

i   X- 

olo  =  .  II:l-7  0IC7         W0  -  T^ 

(12)        *{  *i  =  '  5  */,-  "Sp 

»■  »■  IS 

Thus,    (11)   becomes 


v/here 

^m  ^<  -    \.  1  -745^(^64       '       (£>-   ."2.^54-0  ^.34 

It  can  be  shown  that  the  error  is  of  the  order  of 


cm  o^^[*(e>"*(^-A(ft)-(^)ui!  vWo^ 


3- 

v;here  f  represents  the  integrand  in  (9)«  The  magnitude  of  the 
relative  error  for  Cu+  is  plotted  in  Graph  1 „   It  is  never 
greater  than  .25% 

To  use  (15)  it  is  necessary  to  develop  an  interpolation 

'Lanczos,  C„,  Applied  Analysis,  Prentice  Hall,  N.J.,  195^* 
pgs  4o4-7 
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formula  for  calculating  the  5  numerators  on  the  right  from  the 
available  table  of  values  of  Z*   This  has  been  done  using  a 
Gregory-Newton  second  degree  divided  difference  interpolation 
formula.  A  third  degree  formula  which  would  be  consonant  with 
the  basic  assumption  involved  in  the  Gaussian  5-point  method 
would  have  increased  the  labor  of  programming  significantly 
and  would  have  made  no  difference  in  the  final  result,  to  the 
limits  of  accuracy  retained. 

The  resulting  three  formulas  are  of  the  form 

Appendix  1  gives  the  CDC-l6o4  program  for  calculating  v, 
starting  at  P  =6  Bohr  radii  and  proceeding  inward.  v(  0  )   is 
calculated  by  equation  (42)  for  0    =5.5  Bohr  radii  as 
-.18181818181818.  The  first  v  calculated  by  the  program  is 
that  for  0     =5  Bohr  radii.  The  values  obtained  are  given  in 
Table  1 . 

Equation  (1)  can  now  be  solved.   It  is  first  desireable  to 
write  it  in  terms  of  Bohr  radii  and  Rydberg  energy  units  to  get 
a  dimensionless  equation. 
Let 

(,7)  e-  £ 

where  a  is  a  Bohr  radius.   Then 

(18)  l^_  ^_L     ^g 


TABLE  I 

( 
Z 

:u+ 

Ag+ 
Z 

ir 

4 

ku+ 

r 

ir 

z 

.000 

29.000 

47.00 

79.00 

.002 

47  c  00 

- 

15910. 

.oo4 

46.99 

- 

7458. 

.005 

28.992 

- 

566Q.2 

.006 

46.95 

- 

4987. 

.008 

46.91 

- 

5705. 

.010 

28.946 

- 

2769.7 

46.85 

- 

2919. 

.015 

28.875 

- 

1805.0 

46.65 

- 

1945. 

.020 

28.759 

- 

1524.4 

46.56 

- 

1458. 

77.55 

-  5502. 

.025 

28.615 

- 

1057.5 

46.06 

- 

1129. 

.050 

28.451 

- 

847.07 

45.78 

- 

920.6 

055 

28.280 

- 

711.96 

45.50 

- 

771.6 

.04 

28.106 

- 

611.27 

45.25 

- 

659.8 

75.55 

-  1575. 

.05 

27.765 

- 

471 .60 

44.75 

- 

509.6 

.06 

27.451 

- 

579.58 

44.20 

- 

409.2 

72.86 

-     955.5 

.07 

27.161 

- 

514.57 

45.62 

- 

557.8 

.08 

26.884 

- 

266.50 

h2.96 

- 

284.7 

70.14 

-     655.6 

.09 

26.608 

- 

229.14 

42.26 

- 

245.8 

.10 

26.520 

- 

199.71 

41.52 

- 

211.6 

67.76 

-    485.1 

.11 

4o.7S 

- 

185.7 

.12 

25.689 

- 

156.27 

4o.o4 

- 

164.5 

65.41 

-     571.9 

.14 

24.977 

- 

126.05 

58.71 

- 

152.8 

62.92 

-    295.5 

.16 

24.209 

- 

104.07 

57.55 

- 

109.5 

60.51 

-    24o.5 

.18 

25.415 

- 

87.555 

56.48 
7 

- 

91.81 

57.75 

-     199.4 

TABLE 

I   (Continued) 

z 

IT 

Ag 

IT 

kx 

+ 

r 

Z 

*> ■ 

Z 

IT 

.20 

22.646 

-  74.745 

55.48 

-  11.99 

55.44 

-  168.0 

.22 

21.929 

-  64.625 

54. 48 

-  66.97 

55.4o 

-  145.5 

.24 

21.279 

-  56.452 

55.45 

-  58.07 

51.66 

-  125.4 

.26 

20.699 

-  49.755 

52.51 

-  50.78 

50.09 

-  107.1 

.28 

20.187 

-  44.124 

51.16 

-  44.76 

48.62 

-     95.55 

•  50 

19.756 

-  59.579 

29.98 

-  59.75 

47.15 

-     82.14 

.52 

28.80 

-  55.55 

45.62 

-     72.47 

.54 

27.66 

-  51.95 

44.05 

-    64.21 

.55 

18.791 

-  50.224 

.56 

26.57 

-  28.89 

42.56 

-     57.15 

.58 

25.55 

-  26.24 

40.65 

-     51    08 

.40 

17.964 

-  25.650 

24.61 

-  25.96 

58.94 

-     45.84 

.^5 

17.114 

-  18. 765 

22.65 

-  19.48 

54.75 

-     55.65 

.50 

16.195 

-  15.050 

21.16 

-  16.06 

50.96 

-     28.41 

.55 

15.196 

-  12.189 

20.04 

-  15.57 

27.82 

-     25.12 

.60 

14.155 

-     9.965 

19.15 

-  11.25 

25.55 

-     19.12 

.65 

25.54 

-     16.05 

.7 

12.051 

-    6.851 

17.58 

-     8.124 

21.75 

-     15.57 

.8 

10.119 

-     4.890 

15.94 

-     5.994 

19.25 

-      9.958 

.9 

8.457 

-     5.616 

14.21 

-     4.505 

17.15 

-      7.419 

1.0 

7.110 

-     2.765 

12.50 

-     5.451 

15.19 

-      5.625 

1.1 

6.018 

-     2.175 

10.87 

-     2.695 

15.54 

-      4.554 

1.2 

5.148 

-     1.758 

9.4o 

-     2.15 

11.68 

-      5.591 

1.5 

4.455 

-     1.454 

10.18 

-      2.695 
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TABLE 

I   (Continue 

,d) 

i 

3u+ 

Ag" 

h 

or" 

Au 

Z 

i- 

r 

z 

Z 

nr- 

1.4 

5.896 

- 

1.228 

7.05 

- 

1.47 

8.88 

- 

2.14 

1.5 

7.76 

- 

1.78 

1.6 

5.070 

- 

.9275 

5.51 

- 

1.08 

6.78 

- 

1.48 

1.8 

2.507 

- 

.7587 

4.11 

- 

.855 

5*25 

- 

1.08 

2.0 

2.107 

- 

.6152 

5.26 

- 

.674 

4.14 

- 

.825 

2.2 

1.817 

- 

.5256 

2.65 

- 

.566 

5.54 

- 

.658 

2.4 

1.6o6 

- 

.4517 

2.21 

- 

.490 

2.72 

- 

.547 

2.6 

1.448 

- 

.4155 

1.89 

- 

.454 

2.28 

- 

.458 

2.8 

1.550 

- 

.5755 

1.65 

- 

.591 

1.95 

- 

.411 

5.0 

1.245 

- 

.5451 

1.48 

- 

.558 

1.70 

- 

.568 

5.2 

1.179 

- 

.5201 

1.55 

- 

.551 

1.52 

- 

.555 

5.4 

1.150 

- 

.2990 

1.25 

- 

.505 

1.58 

- 

.509 

5.6 

1.095 

- 

.2808 

1.19 

- 

.285 

1.28 

- 

.287 

5.8 

1.068 

- 

.2651 

1.15 

- 

.267 

1.21 

- 

.269 

4.0 

1.050 

- 

.2512 

1.10 

- 

.252 

1.15 

- 

.255 

4.5 

1.022 

- 

.2226 

1.04 

- 

.225 

1.07 

- 

.225 

5.0 

1.008 

- 

.2001 

1.02 

- 

,200 

1.05 

_ 

.200 

5.5 

1.005 

- 

.1818 

1.01 

- 

.182 

1.01 

- 

,182 

6.0 

1.001 

— 

.1667 

1.00 

_ 

.167 

1.00 

_ 

.167 

so  that  (1)  becomes 

(19)  -p-V4"  —  G=-V)E  =  o 

If  E  is  in  Rydberg  units  we  can  write 

(20)  *'-  —£Z    *  , 

where  £  is  a  dimensionless  number,  since  a  Rydberg  unit  is 
given  by 


(21) 

2-  *^ 

Also,    if 

(22) 

^P 
or  r            r 

then 

v  *w(£) 

(23) 

Inserting  (20)  and  (25)  in  (19), 

But, 
(25)  0.,=  Jll 


r^  €' 


so  (24)  becomes 

(26)  ^i£  1-  I  c  -  £(?)>=  o 


10 


where 


(27) 


4=  a 


IT 


Equation  (26)  can  be  solved  readily  using  a  Runge-Kutta 
4th  order  approximation.   Let 


Then 

and  the   set  of  equations  to  be   solved  is 


gift.  _  J>(2 


(50) 


<&--i>(fV«l* 


C3tps  * 

The  formulas  for  the  solution  are 

where  the  quantities  on  the  right  are  calculated  in  the  following 
order: 


(52) 


Here, 


(55) 


k*  e,-Po 


11 


After  making  the  calculations  (52)  and  inserting  in 
formulas  (51)  we  get 

*  «(^)^  J^^)-#&H} 

As  before,  <Q  (^  — -jis  calculated  using  a  Gregory- 

Hevrton  second  degree  divided  difference  interpolation  formula 
of  the  same  form  as  equation  (16). 

Since  equations  (^4)  give  only  successive  increments  in 
R  and  y  it  is  necessary  to  obtain  starting  values  for  each. 

Assume  that  near  the  origin 

(55)         e.Wo^V^eV ..  . , 

Also,  for  small  £>    , 

(56)  £=   -  i-£ 

r 

where  Z  is  the  complete   charge  on  the  ion.      Differentiating 
(55)   tv/ice  and  inserting  it  and   (56)   in  (26)   gives 


(57) 

Thus , 
(58) 


4-    iiiai0 


^*3^  . 
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For   small  ,   we  want 

09)  ;  ^      S  f 

Inserting  (39)  in  (35) 

(4o)    R(S?)=  b,(Sf)  +  bi  #b,(S(^ 

so  that  b,  -  1  • 

As  starting  values,  let  us  calculate  initial  values  for 
R  and  y  at  0   =.002.   Inserting  the  coefficients  in  (35) 

V,re  can  ignore  £  compared  with  2Z  .   Then  the  starting 
values  for  any  £  are 

Cu+  Ag+  Au+ 

R  .002  .002  .002 

y  .887         .821        .709 

The  CDC-15o4  program  for  the  calculation  is  in  Appendix  2. 
The  magnitude  of  the  initial  energy  was  set  at  .4  in  each  case 
and  increased  in  increments  of  .01  to  1.0  (to  1.6  for  gold). 
The  potentials  for  (°  =.002  were  calculated  from  equation  (36) • 
Potentials  f or  P>  6  were  calculated  from 

(42)  r  s  -L      . 

The  wave  functions  from  this  calculation  are  intended 

merely  for  starting  values  in  equation  (3),  so  no  great  accuracy 

8  1  5 

is  required.  The  error  is  of  the  order  of  h  .  Only  a  sample 

°Kunz,  K.,  op.  cit.,  pg  188 
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curve  of  R  vs  P  is  shown  for  each  metal,  that  for  the  energy 
value  at  which  the  minimum  of  the  curve  of  fe.  vs  ^  was  achieved. 
This  is  given  in  Graph  2.   The  curves  of  c  vs  ^  are  shown  in 
Graph  J. 

To  solve  (5)  it  may  first  be  put  in  dimensionless  form: 


A^ 


v/here 


(44) 


G(f,r)»£  2  :        .,»(?)  WW 


*x    X 


and  r  *-l 

)       "^ 

Here  the  f  .  refer  to  the  wave  functions  for  the  individual 
n,l 

electronic  shells.  These  are  available  (cf  note  4). 

For  Cu+,  the  filled  shells  are  the  1s,  2s,  2p,  ps,  5p>  and 
3d,  corresponding  to  n,l  values  of  10,  20,  21,  JO,  51  >  52.  The 
integral  then  becomes  (using  r  *£.  P  as  an  example) 


no 


p     /0  5l'  ^fT  -J(J     ^l  ' 

and  similarly  for  r  > /9  , 

All  the  integrals  on  the  right  can  be  evaluated  using 
numerical  methods  with  the  wave  functions  cal<  ulated  above  as 
starting  values,  although  the  programming  will  be  tedious. 
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The  integ? ation,  for  all  practical  purposes,  probably  need  not 
proceed  beyond  the  range  0  to  6,   The  right  hand  side  of  (46) 
then  becomes 

where  the  a's  are  the  constants  resulting  from  the  integrations. 
It  is  then  necessary  to  solve  the  set  of  equations 

(48)  ; 

by  a  Runge-Kutta  approximation  using  equations  (31 )  and  (p2) 


a^tatfVfelK^ 


as  before.   The  resulting  wave  functions  are  then  used  to 
re-evaluate  the  integrals  on  the  right  of  (46)  and  so  on.   The 
procedure  is  terminated  when  succeeding  sets  of  wave  functions 
differ  by  less  than  the  required  degree  of  accuracy. 

Unfortunately,  there  was  not  enough  time  available  to  write 
the  CDC-1o04  program. 
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CONCLUSIONS 

The  results  obtained  differ  significantly  from  thoseof  other 
investigators.  For  example,  Fuchs  calculated  the  minimum  energy  for 
copper  as  approximately  .8,  occurring  at  about  2.9  Bohr  radii. 
He  gives  no  estimate  of  his  possible  error,  however. 

Since  the  error  in  finding  each  increment  in  the  method 
used  here  is  of  the  order  of  hr  ,  the  cumulative  error  in  R  to  the 
radius  at  which  the  minimum  occurs  is  in  no  case  greater  than  ,005. 
The  accuracy  of  the  method  cannot  be  called  into  question,  therefore, 

There  is  always  the  possibility  of  a  mistake  in  algebra  or 
programming,  such  as  the  omission  of  a  factor  of   two.   Every  step 
in  the  calculation  has  been  checked  over  four  times,  and  if  there  is 
such  an  error  it  must  be  very  subtle  in  nature. 

It  should  not  be  thought,  however,  that  the  results  of 
previous  investigators  in  this  field  are  necessarily  correct.  They 
were  calculating  by  hand,  and  the  methods  which  they  were  forced  to 
use  were  so  oversimplified  that  reasonable  accuracy  v/as  often 
impossible.   Indeed,  if  the  previous  investigators  (who  all  seem 
to  have  used  the  same  methods)  had  been  in  error  only  .01  units 
for  each  increment  of  R,  their  total  error  could  be  very  large. 
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APPENDICES 


Note:  Figures  in  the  2  column  refer  to  A^  + 
and  in  the  5  column  to  Au  ,  where  different 
from  the  1^°  column,  which  is  for  C\. 
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APPENDIX  I 


ORG  10000 

ST OR    EQU  11000 

MS     EQU  1 2000 

VR1     EQU  11100 

LDA  0  VII! 

3TA  0  VR1 

SFT     ENI  1  00000 

ENI  2  00000 

BEGIN   LDA  0  ALPHA 

FMU  1  RAD+1 

STA  0  STOR 

LDA  0  BETA 

PMU  1  RAD+2 

FAD  0  STOR 

STA  0  STOR 

FDV  0  TWO 

STA  0  STOR+1 

LDA  0  STOR 

FMU  0  STOR 

STA  0  STOR 

LDA  0  BETA 

FMU  1  RAD+1 

STA  0  STOR+2 

LDA  0  ALPHA 

FMU  1  RAD+2 
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FAD  0  STOR+2 

STA  0  STOR+2 

FDV  0  TWO 

STA  0  STOB+2 

LDA  0  STOR+2 

FMU  0  STOR+2 

STA  0  STOR+2 

LDA  1  RAD+1 

FAD  1  RAD+2 

STA  0  STOR+4 

FDV  0  TWO 

STA  0  STOR+5 

LDA  0  STOR+4 

FMU  0  STOR+4 

STA  0  STOR+4 

LDA  1  ZSTA 

FSB  1  ZETA+1 

STA  0  STOR+6 

LDA  1  ZETA+1 

F3  1  ZETA+2 

STA  0  STOR+7 

LDA  1  RAD 

FSB  1  RAD+1 

STA  0  3T0R+8D 

LDA  1  RAD+1 

FSB  1  RAD+2 

STA  0  STOR+9D 
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LDA  1  RAD 

FSB  1  RAD+2 

ST A  0  ST0R+10D 

LDA  0  STOR+6 

FLU  0  3T0R+9D 

3TA  0  ST0R+11D 

LDA  0  STOR+7 

FMU  0  ST0R+8D 

3TA  0  3T0R+12D 

LDA  0  3TCR+11D 

FSB  0  3T0R+12D 

STA  0  STOR+1 1D 

LDA  0  ST0R+8D 

FMU  0  3T0R+9D 

FMU  0  ST0R+10D 

STA  0  3T0R+9D 

LDA  0  3T0R+1 

FSr  1  RAD 

STA  0  STOR+1 JD 

FDV  0  ST0R+3D 

FMU  0  STOR+6 

STA  0  3T0R+UD 

LDA  0  STOR+1 

FSB  1  P"D+1 

FMU  0  STOR+1 pD 

FDV  0  3TCR+10D 

FMU  0  STOR+1 1D 

FAD  0  3TOR+14D 

FAD  1  ZETA 

FDV  0  3T0R 

STA  0  STOR+1 
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LDA  0  STOR+5 

FSB  1  RAD 

STA  0  STOR+1 5D 

FDV  0  3T0R+8D 

FMU  0  STOR+6 

STA  0  ST0R+16D 

LDA  0  STOR+3 

FSB  1  RAD+1 

FMU  0  STOR+1 5D 

FDV  0  5T0R+10D 

FMU  0  3T0R+11D 

FAD  0  ST0R+16D 

FAD  1  ZETA 

FDV  0  STOR+2 

FAD  0  STOR+1 

FLU  0  FIVE 

STA  0  3TOR+5 

STA  0  STOR+5 

LDA  0  3T0R+5 

FSB  1  RAD 

STA  0  3TOR+17D 

FDV  0  3T0R+8D 

FMU  0  3T0R+6 

STA  0  STOR+1 8D 

LDA  0  STOR+5 

FSB  0  RAD+1 


2p 


PMU  0  ST0R+17D 

FDV  0  3T0R+10D 

PMU  0  3T0R+11D 

FAD  0  ST0R+18D 

FAD  1  ZSTA 

fdv  0  STOR+4 

FMU  0  FIGHT 

FAD  0  3T0R+5 

FMU  0  TWO 

FDV  0  NINE 

STA  0  3TOR+5 

LDA  1  RAD+2 

FSB  1  RAD+1 

FMU  0  3TOR+5 

FAD  0  VR1 

STA  0  VR1 

LDQ  1  RAD+2 

STQ  2  ANS 

STA  2  ANS+1 

INI  2  00002 

ENI  0  00000 

I3K  1  55D  ISK  1   59D  ISK  1  49D 


3LJ  0  BEGIN 

ENI  5  00000 

ENI  4  ANS 

OONV   SLJ  4  71000 
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Etfl  0  00000 

001  0  ANS 

006  0  0 

INI  4  00006 

SNA  4  00000 

SAU  0  C0NV+1 

ENI  0  00000 

ISK  5  17D  ISK  5   19D  I3K  5   16D 

SLJ  0  00NV 

ENI  0  00000 

SL3  0  00000 
VIN   DSO    -.18181818181818 
ALPHA  DEC    1 .77^596666 
BETA  DSC    .22540554 
TWO   DEC    2. 
FIVE  DEC    5. 
EIGHT  DEC    8. 
NINE  DEC    9. 
RAD   DEC    6. 

DEC    5.5 

DEC    5. 

DEC    4.5 

DEC    4. 

DEC    5.8 

DSC    5.6 

DEC    5.4 
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DSG 

5.2 

DSC 

5. 

DEC 

2.8 

DEC 

2.6 

DEC 

2.4 

DEC 

2.2 

DEC 

2. 

DEC 

1.8 

DEC 

1.6 

DEC 

1.4 

DS3 

1.5 

DEC 

1.2 

DEC 

1.1 

DEC 

1. 

DSC 

.9 

DEC 

.8 

DEC 

.7 

DEC 

.6 

DEC 

.55 

DSC 

.5 

DSC 

.45 

DEC 

.4 

DEC 

.55 

DEC 

.5 

DSC 

.28 

DEC 

.26 

DEC 

1.4 

DEC 

1.2 

DEC 

1.1 

DSC 

1.0 

DEC 

.9 

DEC 

,8 

DEC 

.7 

DEC 

.6 

DEC 

.55 

DEC 

.5 

DEC 

.45 

DEC 

.4 

DEC 

.58 

DSC 

.56 

DSC 

.54 

DSC 

.52 

DSC 

.5 

DSC 

1.5 

DEC 

1.4 

DSC 

1.5 

DSC 

1.2 

DEC 

1.1 

DEC 

1. 

DEC 

.9 

DEC 

.8 

DEC 

.7 

DEC 

.65 

DEC 

.o 

DSC 

.55 

DEC 

.5 

DEC 

.^5 

DEC 

.4 

DEC 

.58 

DEC 

.56 
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DEO 

.24 

DEC 

.28 

DSC 

.5* 

DEC 

.22 

DSC 

.26 

DEC 

.52 

DEC 

.2 

DEC 

.24 

DEC 

.5 

DEC 

.18 

DEC 

.22 

DEC 

.28 

DEC 

.16 

DEC 

.2 

DEC 

.26 

DEC 

.14 

DEC 

.18 

DEC 

.24 

DEC 

.12 

DEC 

.16 

DEC 

.22 

DEC 

.1 

DEC 

.14 

DEC 

.2 

DEC 

.09 

DEC 

.12 

DEC 

.18 

DEC 

.08 

DEC 

.11 

DEC 

.16 

DEC 

.07 

DEC 

.1 

DSC 

.14 

DEC 

.06 

DEC 

.09 

DEC 

.12 

DEC 

.05 

DEC 

.08 

DEC 

.1 

DEC 

.04 

DEC 

.07 

DEC 

.08 

DEC 

.055 

DEC 

.06 

DEC 

.06 

DEC 

.05 

DEC 

.05 

DEC 

.04 

DEC 

.025 

DEC 

.04 

DEC 

.02 

DEC 

.02 

DEC 

.055 

DEC 

0 

DEC 

.015 

DSC 

.05 

ZETA  DSC 

1. 

DSC 

.01 

DEC 

.025 

DEC 

1.01 

DSC 

.005 

DEC 

.02 

DSC 

1.05 

DSC 

0 

DEC 

.015 

DEC 

1.07 

ZETA 

DSC 

1.001 

DSC 

.01 

DEC 

1.15 

DEC 

1.005 

DEC 

.008 

DEC 

1.21 

DEC 

1.008 

DEC 

.006 

DEC 

1.28 

DEC 

1.022 

DSC 

.oo4 

DEC 

1.58 
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DEO  1.050  DEO   .002  DEC  1.52 

DEC  1.068  DEC   0  DEC  1.70 

DEC  1.095        ZETA  DEC   1  DEC  1.95 

DEC  1.150  DEC   1.01  DEC  2.28 

DEC  1.179  DSC   1.02  DSC  2.72 

DEC  1.245  DEC   1.0 A  DEC  5.54 

DEC  1.550  DEC   1.10  DSC  4.14 

DSC  1.448  DEC   1.15  DEC  5.25 

DEC  1.606  DEC   1.19  DEC  6.78 

DSC  1.817  DSC   1.25  DEC  7.76 

DEC  2.107  DEC   1.55  DEC  8.88 

DEC  2.507  DEO   1.48  DEC  10.18 

DEC  5.070  DSC   1.65  DSC  11.68 

DEC  5.896  DSC   1.89  DSC  15.54 

DSC  4.455  DEC   2.21  DEC  15,19 

DEC  5.148  DEC   2.65  DEC  17.15 

DEC  6.018  DEC   5.26  DEC  19.25 

DSC  7.110  DEC   4.11  DEC  21.75 

DSC  8.467  DEC   5.51  DEC  25.54 

DSC  10.119  DEC   7.05  DEC  25.55 

DSC  12.051  DEC   o.4o  DEC  27.82 

DEC  14.155  DEC  10.87  DEC  50.96 

DEC  15.196  DSC  12.5  DEC  54.75 

DEC  15.195  DEC  14.21  DEC  58.94 

DEC  17.114  DEC  15.94  DEC  40,65 

DEC  17.964  DEC  17.58  DEC  42.56 
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DEC  18.791 

DEO  19.7^6 

DEC  20.187 

DEC  20.699 

DEC  21.279 

DEC  21.929 

DSC  22.645 

DEC  25.415 

DEC  24.209 

DEC  24. 977 

DEC  25.689 

DEC  26.520 

DEC  26.608 

)  26.884 

DEC  27.161 

DSC  27.451 

DEC  27.765 

DSC  28.106 

DEC  28.280 

DEC  28.451 

DEC  28.615 

DSC  28.759 

DSC  28.875 

DEC  28.9^8 

DEC  28.992 

DEC  29.000 

s;:d 


DSC  19.15 

DEC  20.04 

DEC  21.16 

DEC  22.65 

DEC  24.61 

DEC  25.55 

DEC  26.57 

DEC  27.66 

DEC  28.8 

DEC  29.98 

DEC  51.16 

DEC  52.51 

DEC  55.45 

DEC  54.48 

DEC  55.48 

DEC  56.48 

DEC  57.55 

DEC  58.71 

DEC  4o.o4 

DEC  40.78 

DEC  41.52 

DSC  42.26 

DEC  42.96 

DEC  45.62 

DEC  44.20 

DSC  44.75 

DSC  45.25 

29 


DEC  44.05 

DEC  45.62 

DEC  47.15 

DEC  48.62 

DEC  50.09 

DEC  51.66 

DEC  55.4o 

DEC  55.44 

DEC  57.75 

DEC  60.51 

DEC  62.92 

DEC  65.41 

DEC  67.76 

DEC  70.14 

DEC  72.86 

DEC  75.55 

DSC  77.55 

DEC  79. 

e::d 


DEC  45.5 

DEC  45.73 

DEC  46.06 

DEC  45.56 

DEC  46.65 

DEC  46.85 

DEC  46.91 

DEC  46 .  95 

DEC  46.99 

DEC  47 . 

DEC  4-7 . 
END 

Note  1 :  Library  subroutine  DECOF  1  must  be  placed  in  memory 
starting  at  cell  71000.  Output  will  be  on  magnetic 
tape  unit  4. 

"ote  2:  Program  prints  P  and  V .  Two  to  a  line  in  decimal 
floating  point. 
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APP21TDIX     II 


REM 

SCHR  EQ  B 

y  RK  METHOD 

ORG 

20000 

VALUE 

EQU 

22000 

3T0R 

EQU 

21000 

\::z 

EQU 

21600 

ARE 

EQU 

21650 

ENERGY 

EQU 

21700 

LDA 

0 

INEN 

3TA 

C 

EIIERGY 

Ell  I 

00000 

Ml 

0 

00000 

;  in 

ENI 

1 

00000 

ENI 

2 

00000 

LDA 

0 

Ri:; 

3TA 

0 

ARE 

LDA 

0 

yil: 

STA 

0 

WYE 

RETURN 

LDA 

1 

RAD 

FSB 

1 

RAD+1 

STA 

0 

3T0R 

LDA 

1 

RAD 

FSB 

1 

RAD+2 

STA 

0 

ST0R+1 

LDA 

1 

RAD-i1 

FAD 

1 

RAD-)2 
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FSB  1  RAD 

FSB  1  RAD 

STA  0  STOR+2 

LDA  1  POT+1 

FSB  1  POT+2 

Fl-TU  0  STOR+2 

FDV  0  FOUR 

FDY  0  STOR+1 

STA  0  STOR+5 

LDA  1  POT 

FSB  1  POT+1 

FMU  0  STOR+2 

FMU  0  3T0R+2 

FDY  0  FOUR 

FDV  0  STOR 

FDV"  0  STOR+1 

STA  0  STOR+4 

LDA  1  POT 

F3.  0  STOR+4 

FAD  0  3TOR+5 

STA  0  STOR+5 

LDA  1  RAD+1 

PSE  1  RAD 

STA  0  STOR+6 

FMU  0  STOR+6 

STA  0  STOR+7 
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FMU  0  STOR+6 

3TA  0  ST0R+8D 

IDA  1  POT 

FSB  0  ENERGY 

STA  0  STOR+9D 

LDA  0  STOR+5 

FSB  0  ENERGY 

FMU  0  STOR+9D 

EMU  0  ST0R+8D 

J  0  ARE 

FDV  0  FOUR 

STA  0  3T0R+10D 

LDA  0  ENERGY 

I  0  WYE 

FMU  0  STOR+7 

STA  0  STOR+11D 

LDA  0  3TOR+5 

FMU  0  WYE 

EMU  0  STOR+7 

STA  0  ST0R+12D 

LDA  0  STOR+5 

FMU  0  ARE 

FMU  0  STOR+6 

FMU  0  TWO 

STA  0  STOR+1?D 

LDA  0  ENERGY 
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FKU  0  THREE 

PMU  0  ARE 

FKU  0  STOR+5 

STA  0  ST0R+14D 

LDA  1  POT 

FMU  0  ARE 

FMU  0  STOR+6 

STA  0  ST0R+15D 

LDA  0  WYE 

FMU  0  SIX 

FAD  0  ST0R+15D 

FSB  0  ST0R+14D 

FAD  0  ST0R+15D 

FAD  0  ST0R+I2D 

FSI  0  3T0R+11D 

FAD  0  ST0R+10D 

FMU  0  STOR+6 

FDV  0  SIX 

STA  0  STOR+16D 

LDA  1  POT+1 

FMU  0  WYE 

FMU  0  STOR+6 

STA  0  ST0R+18D 

LDA  0  WYE 

FMU  0  STOR+6 

FDV  0  TWO 


yh 


FAD  0  ARE 

STA  0  STCR+19D 

LDA  0  3T0R+5 

FSB  0  ENERGY 

STA  0  3T0R+20D 

LDA  1  POT+1 

FSB  0  SIT2RGY 

FMU  0  ST0R+20D 

FMU  0  ST0R+19D 

FMU  0  3TCR+7 

FDV  0  TWO 

STA  0  3T0R+21D 

LDA  0  ENERGY 

PM  0  ARE 

FMU  0  SIX 

STA  0  3T0R+22D 

LDA  1  POT+1 

FMU  0  ARE 

STA  0  STOR+25D 

LDA  1  POT 

FS]  0  ENERGY 

STA  0  ST0R+24D 

LDA  0  3TOR+5 

FSB  0  ENERGY 

FMU  0  STCR+24D 

FMU  0  ARE 
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FMU  0  STCR+7 

FDV  0  TV/0 

3TA  0   ST0R+25D 

LDA  0  E1TZR"Y 

FMU  0  WYE 

FMU  0  STOR+6 

FMU  0  THREE 

STA  0  ST0R+26D 

LDA  0  STOR+5 

fmu  c  wye 

FMU  0   3T0R+5 

FMU  0  TWO 

STA  0  ST0R+27D 

LDA  0   STOR+5 

FMU  0  ARE 

FMU  0  FOUR 

STA  0  ST0R+28D 

LDA  1    POT 

FMU  0  ARE 

FSB  0  3T0R+22D 

F  lD  0  3TOR+28D 

FAD  0    3TOR+27D 
0   3T0R+26D 

FAD  0   ST0R+25D 

FAD  0   SV0R+2^D 

FAD  0   3T0R+21D 

FAD  0   ST0R+18D 


PMU  0  3T0R+6 

FDV  0  SIX 

FAD  0  WYE 

STA  0  WYE 

3TA  2  VALUE+5 

LDA  1  RAD+1 

STA  2  VALUE 

LDA  0  3T0R+1  'D 

FAD  0  ARE 

STA  0  ARE 

STA  2  VALUE+1 

FDV  1  RAD+1 

STA  2  VALUE+2 

INI  2  00004 

si:i  0  00000 

I3K  1  5oD  ISK  1  591  ISK  1  52D 

3LJ  0  RETURN 

ENI  4  00000 

ENI  5  VALUE 

CONV  3LJ  4  71000 

ENI  0  COOOO 

001  0  VALUE 

004  0  0 

INI  5  oooo4 

ENA  5  00000 

SAU  0  OONV+1 
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ENI  0  00000 

ISK  4  56D  ISK  4  590  ISK  4  52D 

3LJ  0  CONV 

LDA  0  ENERGY 

FAD  0  DELEN 

STA  0  ENERGY 

ENI  0  00000 

ISK  y   00075  ISK  5  00170 

SLJ  0  BEGIN 

ENI  0  00000 

SL3  0  00000 

INEN   DEC  -.4 

TWO    DEC  2. 

THREE  DEO  3. 

FOUR   DEC  4. 

SIX    DEO  6, 

DELEN  DEC  -.01 

RIN    DEC  .002 

YIN    DEC  .887  DEC  .821  DEC  .709 

RAD    DSC  .002 

DEC  .005  DEC   .004  DEC  .02 

DEC  .01  DSC  .006  DEC  .04 

DEC  .015  DSC   .008  DSC   .06 

DEC  .02  DEC   .01  DEC  .08 

DSC  .025  DSC   .015  DSC   .1 

DEC  .05  DEC   .02  DEC  .12 
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DEC 

.055 

DEC 

eo4 

DEO 

.05 

DEC 

.06 

DEC 

.07 

DEC 

.08 

DEC 

.09 

DEC 

.1 

DEC 

.12 

DEC 

.14 

DEC 

.16 

DEC 

.18 

DSC 

.2 

DSC 

.22 

DSC 

.24 

DSC 

.26 

DSC 

.28 

DEC 

.5 

DEC 

.55 

DEC 

.4 

DSC 

.45 

DEC 

.5 

DEC 

.55 

DEC 

.6 

DEC 

.7 

DSC 

.8 

DEC 

.025 

DSC 

.05 

DEC 

.055 

DEC 

.04 

DEC 

.05 

DEC 

.06 

DEC 

.07 

DEC 

.08 

DSC 

.09 

DEC 

.1 

DSC 

.11 

DEC 

.'2 

DEC 

.14 

DSC 

.16 

DSC 

.18 

DEC 

.2 

DEC 

.22 

DEC 

.24 

DEC 

.26 

DEC 

.28 

DEC 

.5 

DSC 

.52 

DEC 

.54 

DSC 

.56 

DEC 

.58 

DSC 

.4 

DEC 

.14 

DEC 

.16 

DEC 

0I8 

DEC 

.2 

DEC 

.22 

DEC 

,24 

DEC 

.26 

DEC 

.28 

DEC 

.5 

DEC 

o52 

DEC 

.54 

DEC 

.56 

DSC 

.58 

DEC 

.4 

DEC 

.45 

DEC 

.5 

DSC 

o55 

DEC 

.6 

DEC 

.65 

DEC 

.7 

DEC 

.8 

DEC 

.9 

DEC 

1. 

DEC 

1.1 

DSC 

1.2 

DEC 

1.3 
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DEC 

.9 

DSC 

.45 

DEC 

1. 

DSC 

.5 

DEC 

1.1 

DEC 

.55 

DEC 

1.2 

DSC 

.6 

DEC 

1.5 

DSC 

.7 

DSC 

1.4 

DEC 

.8 

DEC 

1.6 

DEC 

.9 

DSC 

1.8 

DEC 

1. 

DSC 

2. 

DEC 

1.1 

DSC 

2.2 

DEC 

1.2 

DSC 

2.4 

DEC 

1.4 

DEC 

2.6 

DEC 

1.6 

DEC 

2.8 

DSC 

1.8 

DEC 

3. 

DEC 

2. 

DEC 

3.2 

DEC 

2.2 

DEC 

3.4 

DEC 

2.4 

DEC 

5.6 

DEC 

2.6 

DEC 

3.8 

DSC 

2.8 

DEC 

4. 

DEC 

3. 

DSC 

4.5 

DEC 

3.2 

DSC 

5. 

DSC 

5.4 

DSC 

5.5 

DEC 

5.6 

DEC 

. 

DSC 

3.8 

DSC 

6.5 

DSC 

0 

DEC 

7. 

DEC 

4.5 

POT  DEC 

-29000. 

DEC 

5. 

POT 


DEC 

1.4 

DEC 

1.5 

DSC 

1.6 

DSC 

1.8 

DEC 

2. 

DEC 

2.2 

DSC 

2.4 

DSC 

2.6 

DSC 

2.8 

DEC 

3. 

DEC 

3.2 

DEC 

3o4 

DEC 

3«6 

DSC 

5.8 

DEC 

4. 

DEC 

4*5 

DSC 

5. 

DEC 

5.5 

DEC 

6 

DEC 

6.5 

DEC 

7. 

DEC 

-79000. 

DSC 

-  7004o 

DSC 

-  5146. 

DEC 

-  1907. 

DEC 

-  1311. 

ho 


DEO 

-11556\ 

DEC 

5.5 

DEC 

- 

5559.4 

DEC 

0. 

DEC 

- 

5610. 

DEC 

0.5 

DEC 

- 

2648.8 

DSC 

7. 

DEC 

- 

2074.6 

POT  DEC 

-27820. 

DEC 

- 

1694.1 

DEC 

-14880. 

DEC 

- 

1425.9 

DEC 

-  9974. 

D2C 

- 

1222.5 

DEC 

7410. 

DEC 

- 

945.2 

DSC 

-  5858. 

DEC 

- 

759.16 

DSC 

-  5890. 

DEC 

- 

629.14 

DSC 

-  2876. 

DEC 

- 

552.6 

DSC 

-  2258. 

DEC 

- 

458.28 

DSC 

-  1841. 

DSC 

- 

599.42 

DEC 

-  1545. 

DEC 

- 

512.54 

DSC 

-  1520. 

DEC 

- 

252.1 

DSC 

-  1019. 

DEC 

- 

208.14 

DSC 

-  818.4 

DEC 

- 

175.07 

DSC 

-  675.6 

DEC 

- 

149.49 

D2C 

-  569.4 

DSC 

- 

129.25 

DSC 

-  487.6 

DEC 

- 

112.9 

DEC 

-  425.2 

DEC 

- 

99.456 

DEC 

-  571.4 

DEC 

- 

88.246 

DEC 

-  529. 

DEC 

- 

78.758 

DSC 

-  265.5 

DEC 

- 

50.448 

DEC 

-  219. 

DEC 

- 

47.5 

DSC 

-  185.5 

DSC  -966.2 

DEC  -745=8 

DEC  -590.6 

DEC  -480.6 

DEC  -598.8 

DEC  -556. 

DEC  -286.6 

DSC  -245.8 

DEC  -214.2 

DEC  -187.1 

DEC  -154,5 

DEC  -144.9 

DSC  -128.4 

DSC  -114.5 

DEC  -102,2 

DEC  -  91.68 

DSC  -  71.5 

DEC  -  56,32 

DEC  -  46.24 

DEC  -  58o 24 

DEC  -  52.06 

DEC  -  27.14 

DEC  -  19.88 

DSC  -  14084 

DEC  -  11.25 

DSC  -  8.66c 


41 


DSC 

-  57.526 

DEC 

-  50.100 

DEC 

-  24.578 

DEC 

-  19.926 

DEC 

-  15.702 

DEC 

-  9.78 

DEC 

-  7.252 

DEC 

-  5.526 

DSC 

-  4.55 

DSC 

-  5.516 

DEC 

-  2.908 

DSC 

-  2.456 

DSC 

-  1.855 

DSC 

-  1.477 

DEC 

-  1.226 

DEC 

~  1.051 

DSC 

-   .9254 

DSC 

.8266 

DEC 

.7510 

DSC 

.6902 

DSC 

.6402 

DSC 

-  -59B 

DEC 

.5^16 

DCC 

-   .5502 

DSC 

.5024 

DSC 

.4452 

DEC 

- 

156. 

DEC 

-  6.782 

DEC 

- 

155.9 

DSC 

-  5.59 

DEC 

- 

116.1 

DSC 

-  4856 

DEC 

- 

101.6 

DEC 

-  5.56 

DSC 

- 

89.52 

DEC 

-  20 96 

DEC 

- 

79.5 

DEC 

-  2,16 

DSC 

- 

71.06 

DEC 

-  1.65 

DEC 

- 

65.9 

DEC 

-  1.52 

DEC 

- 

57.78 

DSC 

-  1.09 

DEC 

- 

52.48 

DSC 

-  .956 

DEC 

- 

47.92 

DEC 

-  .822 

DEC 

- 

58.96 

DSC 

-  .756 

DSC 

- 

52.12 

DSC 

-   .67 

DEC 

- 

26.74 

DEC 

-   .618 

DSC 

- 

22. 46 

DEC 

-   .574 

DEC 

- 

16.25 

DEC 

-  .558 

DSC 

- 

11.99 

DEC 

-  .506 

DEC 

- 

9.01 

DSC 

-  a446 

DEC 

- 

6,902 

DSC 

-  ,4 

DSC 

- 

5.59 

DEC 

-  .564 

DSC 

- 

4.5 

DEC 

-  .55^ 

DSC 

- 

2.94 

DEC 

-  .508 

DEC 

- 

2.16 

DEC 

-  c286 

DSC 

- 

1 .666 

END 

DSC 

- 

1.548 

DEC 

- 

1.152 
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DEC 

-   .4002 

DEO 

-    .5656 

DEC 

-   .5554 

DEC 

-    .5077 

DEC 

-   .2357 

si;d 

DSC 

-  .98 

DEC 

-  .868 

DEC 

-  .782 

DSC 

-  .716 

DEO 

-  .662 

DEC 

f  4 

DEC 

"      .57 

DEC 

-  c554 

DEO 

-  .504 

DEC 

-  Ah6 

DEC 

-     A 

DEC 

-  .564 

DEC 

-  .554 

DSC 

-  .508 

DEC 

-  .286 

END 

Note  1  :  i'ote  1  to  Appendix  I  also  applies  here 

Note  2:  Program  prints  O     |?  K_  ,x^A   j    in  decimal  floating  point 
along  each  line. 
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